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Method Of Resolving Analytes In a Fluid 

STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER 
FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT 

5 The United States Government may have certain rights to the present 

invention pursuant to Grant No. DAAG55-97-1-0187 awarded by the U.S. Army and 

Grant No. DAAK60-97-K-9503 awarded by the U.S. Navy. 

RELATED CASES 

10 The present patent application claims priority from U.S. Provisional 

Application No. 60/086,876, filed May 27, 1998 and from U.S. Provisional Application 
No. 60/088,804, filed June 9, 1998. 

BACKGROUND OF THE INVENTION 

15 The present invention is related to the electrochemical sensing of analytes 

in a medium or a vacuum, such as odors, smells, vapors and the like in air, and in 
particular, to methods of determining or resolving such analytes once detected. 

There has been considerable interest in developing devices for the general 
detection of analytes in a fluid, vacuum, air, or other medium. A specific instance of an 

20 analyte detector is a device for sensing smells or odors (i.e., analytes in air). It is well 
recognized that some animals, such as dogs, have a much keener sense of smell than 
human beings. Because of the acuity of their noses, dogs have been utilized for many 
tasks including, for example, the detection of bombs, mines, drugs, poison gases, and 
illegal contraband; dogs also aid in the search and rescue of humans. Devices for sensing 

25 smell would be useful for the applications where animals are traditionally used, as well as 
for a multitude of uses where animals are impractical or inappropriate. 

Moreover, devices for the general detection of analytes have potentially 
many more applications than devices for detecting only smells and odors. For example, 
the uses for devices for analyte detection include the detection of chemical leaks, quality 

30 control in food processing, medical diagnosis and testing, fabrication and manufacture of 
commercial and industrial goods, pharmaceutical production, testing or evaluating any 
odorant or analyte in any medium (e.g., fuel, oil, wine, solvents), and many other 
applications. An instrument for analyte detection would be highly desirable in many 
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industries and applications, such as the chemical, petrochemical, food, fragrance, medical, 
automotive, military, environmental, health and safety, and air quality monitoring sectors. 
Therefore, it is desirable to develop techniques and devices for the detection of analytes. 

One approach for sensing analytes is the use of an array of sensors on a 

5 substrate. Each sensor has a chemical compound which has a particular electrical 
characteristic response to exposure to an analyte, such as an odor. The sensors are 
connected to electrical circuits by which the responses of the chemical compounds to the 
analyte are retrieved. See, for example, U.S. Appln. No. 09/130,775, entitled "SENSOR 
ARRAYS ON AN IC FOR ANALYTE DETECTION," filed August 7,1998 and assigned 

10 to the present assignee, and which is incorporated by reference herein. Each site responds 
to a broad class of stimuli, with the collective response of the many different members of 
the sensor array providing a fingerprint of an analyte of interest. The number of data 
points generated for each analyte is large. Even with modern high-speed computers, the 
amount of data generated may significantly slow operations to analyze the sensed 

1 5 analytes under real world constraints, in particular. Therefore, there is a need for 

techniques and systems which are computationally efficient in determining or resolving 
an analyte upon detection from the large amount of data derived from the sensor array. 

The present invention is directed toward such a technique and system for 
quickly and efficiently resolving analytes. A quantitative metric is generated for the 

20 effective resolution of different analytes by the sensor array. The present invention 
further provides a powerful tool in the analysis and construction of olfactory sensor 
arrays. 

SUMMARY OF THE INVENTION 
25 The present invention provides for a method for distinguishing different 

odors by the steps of disposing a plurality of d sensors in an array, each sensor having 
different electrical responses to different orders; exposing the sensors to first and second 
odors; generating first and second sets of data points from each of the sensors, each set 
corresponding to the first or second odor, each data point being represented by a vector in 
30 a d-dimensional space; determining an axis in the d-dimensional space, the axis having 
the property that projections of the data points onto the axis in the d-dimensional space 
have optimal separation; and resolving the first odor from the second odor by the 
separation. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 illustrates an array of sensors connected to a computer unit, in 
accordance with one aspect of the present invention; 

Fig. 2 illustrates a typical electrical response of a sensor to the presence of 

an analyte; 

Fig. 3 is a two-dimensional representation of the resolution factor obtained 
by the projection of the data points of two analytes upon optimal projection vector, 
according to one aspect of the present invention; 

Fig. 4A is a table of pairwise resolution factors obtained from an array of 
fourteen carbon black polymer-composite detectors for nineteen solvent analytes; Fig. 4B 
is a table of pairwise resolution factors obtained from an array of eight tin oxide detectors 
for the nineteen solvent analytes; Fig 4C is a table of pairwise resolution factors obtained 
from an array of twelve bulk conducting organic polymer detectors for the nineteen 
solvent analytes; 

Fig. 5 is a table of normalized pairwise resolution factors obtained from 
the array of fourteen carbon black-polymer composite detectors; 

Fig. 6 is a table of principal components of the fourteen element carbon 
black-polymer composite detector response data; 

Fig. 7 A is a plot of the average ability to resolve all analyte pairs, rf(n) , 
versus the number n of sensors in the array; Fig. 7B is a plot of the average ability to 
resolve all analyte pairs, rf(n) , for the carbon black-polymer composite detectors versus 
the number n of sensors in the array using normalized data; Fig. 7C is a plot of the 
maximum average resolution factor, rf max(n), the specific set of specified number n of 
detectors which have the largest ability on average to resolve all solvent vapor analyte 
pairs versus the number of detectors n; Fig. 7D is a plot of the maximum average 
resolution factor, rfmax(n), for the carbon black-polymer composite detectors versus the 
number of detectors n using normalized data; 

Fig. 8 A is a plot of the smallest pairwise resolution factor, rf WO rsi» versus 
the number n of sensors in the array; Fig. 8B is a plot of the smallest pairwise resolution 
factor, rfwom, versus the number n of carbon black-polymer composite sensors in the 
array; 

Fig. 9 is a plot of the largest solution factor, rf bcs t(n), to resolve ethyl 
acetate from tetrahydrofiiran at fixed concentrations for all combinations of detectors 
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from the number n of the detectors in the array versus the number n of the detectors in the 
array; and 

Fig. 10 is a table of the results of resolving on average and resolving the 
worst-resolved analyte pair using combined arrays of carbon black-polymer composite, 
5 tin oxide, and conducting organic polymer sensors. 

DESCRIPTION OF THE SPECIFIC EMBODIMENTS 
In accordance with the present invention, an array 10 of sensors is 
disposed on a substrate, such as shown in Fig. 1 . Each sensor is responsive to different 

1 0 odor analytes in a particular manner so that the response of all the sensors creates a 
"fingerprint" which can be used to identify the analyte. These sensors may include 
surface acoustic wave (SAW) devices, tin oxide detectors, conducting organic polymers, 
dye-impregnated polymer films on fiber optic detectors, polymer-coated micromirrors, 
quartz crystal microbalances (QCMs), electrochemical gas detectors, chemically sensitive 

1 5 field-effect transistors, carbon black-polymer composite chemiresistors, micro-electro- 
mechanical system (MEMs) devices, and micro-opto-electro-mechanical system 
(MOEMs) devices. If the devices are compatible, the sensors might be placed on a single 
substrate, such as that of semiconductor material, for the benefits of miniaturization and 
convenience. Each of the sensors is connected by wires through an array interface 12 to a 

20 computer unit 1 which receives the electrical responses from the sensors upon exposure to 
an odor analyte. The computer unit is programmed to save and process the data in 
accordance with the data processing steps set forth below. The results may be displayed 
on a screen 5 of the computer unit's monitor and/or saved in memory, such as represented 
by a computer disk 15. 

25 With d sensors 1 1 in the array 10, each sampling of the array creates d 

values, i.e., each sampling data point can be considered a d-dimensional vector. If there 
are eight sensors, each data point has eight values, one from each sensor. Hence it can be 
readily understood that the amount of data generated by the array 10 can be quite large. 
For example, if there are 100 sensors in the array 10, then each data point has 100 

30 associated values. Furthermore, there may be several samples taken for each analyte 
thereby magnifying the computational problem of analysis. The present invention 
reduces the problem of d-dimensional data point vectors to the comparison of values in 
one dimension. This is performed by the application of the Fisher linear discriminant to 
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the data. In this manner the resolution of different odor analytes can be made efficiently 
and effectively with a quantitative metric. 

In the description below, specific types of sensors or detectors (tin oxide, 
conducting organic polymer and carbon black-polymer composite detectors) are 

5 described with respect to the detection of specific analytes (some nineteen different 
solvent vapors including 1 ,2-dimethoxyethane, acetone, acetonitrile, anisole, benzene, 
butylamine, chloroform, cyclohexane, dichloromethane, ethanol, ethyl acetate, 
isopropanol, methanol, n-heptane s n-pentane, tetrahydrofuran, toluene, triethylamine and 
a,a,a-trifluorotoluene). It should be understood that the present invention described 

1 0 herein is applicable to other sensors, including those of different sensing modalities, and 
analytes. 

The present invention might best be understood by assuming a typical 
experiment where it is desired to obtain a quantitative measure of the efficiency of an 
array of 14 detectors in distinguishing between two solvents. The array is exposed to 

1 5 each solvent for 12 cycles. Each cycle consists of 300 seconds of pure carrier gas over 
the detector array, then 240 seconds of gas and the solvent, then 300 seconds of pure 
carrier gas again. This procedure ensures steady state conditions before, during, and 
after, exposure to the solvent mixture. The resistance of each detector, j, is measured at 3 
second intervals throughout each exposure, i. The experiment data base for this typical 

20 experiment thus consists of 2* 1 2*840/3* 14 = 94,080 resistance measurements. 

The first step in processing these data finds the maximum relative 
resistance change during each exposure for each solvent and each detector. Of course, it 
should be understood that the present invention is also applicable to other electrical 
parameters, such as voltage, current, capacitance, or inductance, which might vary in 

25 response to the exposure to analytes. A typical detector response during an exposure will 
have the general form, such as illustrated in Fig. 2. The baseline resistance, Rb, is the 
stable value attained by the detector just prior to the introduction of the solvent into the 
carrier gas stream, that is, the 100 th measurement (300 seconds sampled at 3 second 
intervals). The maximum resistance value attained during the gas and solvent interval, 

30 R m ax, is readily obtained between measurements 101 and 540 (typically towards the end 
of the 240 second gas and solvent interval). The maximum relative resistance change for 
the cycle is then: 

(Rmax - Rb)/Rb- 
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The experimental database has now been reduced to 2*2* 14 = 336 measurements. These 
can be regarded as a set of twelve 14-dimensional points, ay (i = 1 to 12 and j = 1 to 14), 
for each solvent. That is: 
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A quantitative measure of the array efficiency in distinguishing between the solvents is 
obtained using the optimum resolution factor (defined by Fisher's linear discriminant): 



The resolution factor, rf, is from an optimal projection vector, w, where 
w = Sw'^ma- m b ) 

1 5 where m a and mb are the respective 14-dimensional means of each of data sets of the two 
solvents, that is: 
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The within-class scatter matrix, S w , is the sum of the matrices S a arid S b where 

S = E(xij-m)(Xij-m) T 
Each S is the sum over the 12 exposures of 14*14 - dimensional terms, such as: 



{a u -m x a n -m 2 ••• a w -m u ) 



10 Then 

S w = S a + S b ; and 
w = Sw" 1 (m a -m b ) 

d w is the distance along the vector w between the projected means and is computed as: 
d w = [w * (m a -m b )(m a -m b ) T ]^ 

15 

Also, the denominator for //is: 

w T Sw = *JcT 2 a -+al^ 

to obtain the resolution factor, rf, between two solvents. 

20 Stated differently, the Fisher linear discriminant searches for the projection 

vector, w, in the detector space which maximizes the pairwise resolution factor, i.e., rf 9 
for each set of analytes, and reports the value of rf along this optimal linear discriminant 
vector. Fig. 3 is a graphical two-dimensional representation of the data points of analytes 
a and b, and their projected means along the optimal discriminant vector w. The rf value 

25 is an inherent property of the data set and does not depend on whether principal 

component space or original detector space is used to analyze the response data. This 
resolution factor is basically a multi-dimensional analogue to the separation factors used 
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to quantify the resolving power of a column in gas chromatography, and thus the rf value 
serves as a quantitative indication of how distinct two patterns are from each other, 
considering both the signals and the distribution of responses upon exposure to the 
analytes that comprise the solvent pair of concern. For example, assuming a Gaussian 

5 distribution relative to the mean value of the data points that are obtained from the 
responses of the array to any given analyte, the probabilities of correctly identifying an 
analyte as a or b from a single presentation when a and b are separated with resolution 
factors of 1 .0, 2.0 or 3.0 are approximately 76%, 92% and 98% respectively. 

The efficacy of the technique is demonstrated in the empirical data 

1 0 obtained with arrays having different types of detectors, i.e., tin oxide, conducting organic 
polymer and carbon black -polymer composite detectors. In the experiments, the array of 
tin oxide detectors had eight different detectors; the array of conducting organic polymer 
detectors had twelve different detectors; and the array of carbon block-polymer composite 
detectors had fourteen different detectors. Nineteen different solvent vapors (1,2- 

15 dimethoxyethane, acetone, acetonitrile, anisole, benzene, butylamine, chloroform, 

cyclohexane, dichloromethane, ethanol, ethyl acetate, isopropanol, methanol, n-heptane, 
n-pentane, tetrahydrofuran, toluene, triethylamine and a,a,a-trifluorotoluene) were used 

as the analytes. . 

Figs. 4A-4C display the pairwise resolution factors, rf for the three 

20 different types of detector arrays. The table in Fig. 4A presents the resolution factors 
obtained from an array of the fourteen carbon black-polymer composite detectors for all 
171 pairs of the nineteen vapors. The tables in Figs. 4B and 4C display similar data for 
the eight tin oxide detectors and for the twelve available bulk conducting organic polymer 
detectors which had all been exposed to the nineteen test solvent vapors. 

25 These resolution factors allow assessment of the performance of the 

different arrays, or of the performance of subsets of the arrays, in various sensing tasks. 
Two criteria are chosen as measures of array performance: (1) the average ability to 
resolve all analyte pairs, rf , and (2) the ability of the array to resolve the worst-resolved 
analyte pair, rfworst. The rf 's and standard deviations across all analyte pairs for the frill 

30 carbon black-polymer composite, tin oxide, and conducting polymer arrays (Figs. 4A-4C) 
evaluated were 145 ± 93, 27 ± 23 and 18 ± 1 6 respectively, while the rfworst values were 
23, 4.0 and 1.4, respectively. Based on this analysis, the carbon black-polymer composite 
array yielded the largest mean statistical separation of the response patterns produced by 
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exposure of this particular collection of available detectors to this particular collection of 
test solvent vapors exposed at fixed concentrations. Using the same experimental 
conditions, the carbon black composite detector array also yielded the largest separation 
between the worst-resolved pair of analytes. 

5 The present invention also provides a means of determining the resolving 

power of the carbon black-polymer composite detectors at variable analyte 
concentrations. The analysis described above serves as one quantitative measure of the 
classification ability of different detector arrays. However, this metric also includes any 
separation between data clusters that arises from the differences in the amplitudes of the 

10 array response to the vapors in the test analyte set. Two solvents that produce nearly 
identical patterns, but with very different absolute detector response amplitudes, would 
produce two well-separated clusters. The large pairwise rf value for such a solvent pair 
produced by the data analysis method described above is relevant for tasks in which the 
concentration of an analyte is known, or when only detection of changes over time in an 

1 5 otherwise relatively constant analyte composition is of interest. However, such a data 
analysis procedure is clearly not appropriate for assessing array classification 
performance in applications where the concentration of the analyte is not known, or not 
known to be fixed, in advance. 

A measure of array performance under such conditions requires that the 

20 fundamental differences in patterns produced by the various vapors be determined 

without considering differences in the mean magnitude of the detector responses. Carbon 
black-polymer composite detectors containing > 20% by weight carbon black have been 
shown to respond linearly to the concentration of several test analytes over at least an 
order of magnitude of analyte concentration. Hence, a measure of the concentration- 

25 independent ability of the array to classify the various solvents in the test can be obtained 
for the carbon black-polymer composite detectors by normalizing the array responses. 
This procedure is equivalent to adding an extra degree of freedom, analyte concentration, 
to the classification task. 

The normalized response, Aty, of the / h detector in an n detector array to 

30 the of k exposures of a specific analyte is defined as: 
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y.» 

where the summation is over the response of all n detectors to all k exposures of the 
specific analyte. For the carbon black-polymer composite vapor detectors, a comparison 
5 of the 171 pairwise resolution factors obtained from normalized responses (the table in 
Fig. 5) to those obtained only from raw response data (the table in Fig. 4A) shows that 
allowing the vapor concentration to be a floating variable slightly reduces, but does not 
remove, the array's ability to resolve the various test analytes at the concentrations used in 
this study. Thus, the rf decreased from 145 ± 93 to 102 ± 50 when the data were 

10 normalized. This analysis can be used for other detectors which exhibit responses that are 
linear with changes in analyte concentration. 

The empirical data also allows comparison of the resolving power of 
different detector arrays with the same number of detectors at known analyte 
concentrations. This allows comparison between the performance of different types of 

1 5 arrays without possible biases due to differing numbers of detectors in each type of 
system. Thus, the data obtained from the specific carbon black-polymer composite 
detector array having the median rf of all possible combinations of twelve of the fourteen 
carbon black-polymer composite detectors studied are compared to the data produced by 
the available twelve element conducting polymer detector array. Similarly, arrays of 

20 carbon black-polymer composite, bulk organic conducting polymer, and tin oxide 
detectors are compared at a common array size of eight detector elements. In these 
comparisons, unnormalized data are used because the lack of linearity with varying 
analyte concentration of some of the detector modalities prevent meaningful 
normalization of the detector response data in those instances. 

25 The carbon black-polymer composite detector array having the median rf 

of all possible combinations of twelve of the fourteen carbon black-polymer composite 
detectors consist of polymers #2 and 4-14 ( refer to the table in Figure 6). This particular 
array displayed rf = 1 17± 79 for the unnormalized response patterns to the set of test 
analytes at the fixed concentrations. For comparison, the available twelve element bulk 

30 organic conducting polymer array produce rf = 27 ± 23 for this same task. 

The three detector arrays are also compared in resolving power using eight 
element arrays. The rf value using unnormalized data for the eight tin oxide detectors 
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was 27 ± 23, while that for the median-performing eight element carbon black-polymer 
composite array (detectors #3, 6-10, 12 and 13 from the table in Fig. 6) is 62 ± 32 and that 
of the median-performing eight element bulk conducting organic polymer array (detectors 
#1-4, 6, 7, 10 and 1 1 as labeled by the manufacturer Neotronics, Inc.) is 13 ± 10. The 

5 trends in rf values observed for twelve element arrays is thus also present at smaller 
array sizes, for the set of tasks and sets of detectors evaluated. 

Rather than the average resolving power rf 5 the resolution between the 
two poorest-resolved solvent pairs may be used to measure the performance of the arrays. 
To obtain this type of comparison, the array having the median rfworst of all possible 

10 combinations of twelve of the fourteen carbon black-polymer composite detectors was 
compared to the rfworst value of the available twelve detector conducting polymer array. 
The selected twelve element carbon black-polymer array contained detectors #1-3,5 and 
7-14 and displayed rfworst = 20 (for anisole vs. butylamine), whereas the twelve element 
conducting polymer array yielded rfworst "1-4 (for n-heptane vs. a,a,a- 

15 trifluorotoluene). Using a set of only eight detectors allowed a comparison between all 
three types of detector modalities. Using this metric, the median-performing eight 
element carbon black-polymer composite array (detectors #2, 4, 5 and 8-12) had an 
rfworst = 1 1 (anisole vs. butylamine) and the eight element tin oxide array had an rfworst 
of 3.8 (for butylamine vs. triethylamine), while the median-performing eight element bulk 

20 conducting organic polymer array (detectors #1, 2, 4-6, 9, 10 and 12) had an rfworst - 1.0 
(n-heptane vs. ct,a,a-trifluorotoluene). 

The data can also used to determine whether the combination of different 
sensor types realize significant improvement in the resolving power of the test solvents. 
The results of the combined arrays in resolving the nineteen analytes on average, and in 

25 resolving the worst-resolved pair, are shown in the table of Fig. 10. All pairwise 

combinations of complete (fourteen element carbon black-polymer composite, twelve 
element conducting organic polymer, and eight element tin oxide) arrays of different 
sensor type, are investigated, along with the performance of the full combination of all 
three arrays of different sensor types. To analyze the data in a consistent fashion for 

30 different arrays and vapors, autoscaled response data are used to formulate the resolution 
factors. The autoscaled response of the jth sensor to the ith exposure, Ay, is: 
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where |jj and nj represent the mean and standard deviation, respectively, of the maximum 
relative differential resistance response of the jth sensor to the entire group of analytes. 

A slight improvement in rf from 9.7 to 9.8 is observed when the twelve 
5 available conducting organic polymer sensors are added to the fourteen element carbon 
black-polymer composite sensor array. The rf won! t increases more significantly, from 0.91 
to 1.8, when these sensors are added. All other combinations of the full arrays , including 
the addition of the eight tin oxide sensors to the twenty-six sensor carbon black-polymer 
composite and conducing organic polymer sensors, underperform the twenty-six sensor 

10 array's performance for this set of test solvents using autoscaled sensor response data. 
An explanation of the observed performance improvement by adding the bulk organic 
conducting polymer sensors to the carbon black-polymer composite sensors is that the 
latter set is weakest at resolving some polar analytes, while the more polar bulk organic 
conducting polymer sensors are strong in resolving polar analytes. The combination of 

15 the two sensor types, in this case, improves the balance in chemical properties of the 
polymers in the sensor array and thus slightly improves the ability of the array to resolve 
this particular set of analytes. 

The data can also be used to evaluate array performance as a function of 
the number of detectors in the array to determine how many detectors are required for 

20 various tasks. Both rf and rf W orst are evaluated for all possible combinations of one 
through fourteen detectors in the carbon black-polymer composite detector array, one 
through eight detectors in the tin oxide array, and one through twelve detectors in the bulk 
organic conducting polymer array. One goal would be to construct the array that has the 
best average resolution between any pair of vapors that it might encounter. This criterion 

25 is closely related to the goal of constructing the array with the most chemical diversity, so 
that it best separates any set of vapors that is likely to be in the environment. To assess 
the diversity of the detector elements in a given detector modality, for all possible 
combinations of every n-element array of each detector type, the average pairwise 
resolution factor, rf(n) , is computed for all 1 71 different solvent pairs of the nineteen 

30 solvents. The value of rf(n) for each allowed combination is tabulated, and the average 
value of rf(n) for all possible combinations having the specified number of detectors is 
plotted vs. the number of detectors in the array. The computed quantity is thus a 
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statistically-based measure of the average pairwise resolution of the solvents in the 
nineteen vapor test set that would be obtained from the average, unbiased and unsorted 
data arising from an n-element detector array of a given detector modality. 

As displayed in Fig. 7A, these data show that, for all three detector arrays, 
5 the average resolving power of the array increased as the numbers of detectors in the 
array increased. Thus, for situations in which the sensing task is not known in advance, 
such as in quality control applications or in complex environments with varying 
backgrounds, the average performance of these tested detector arrays improved as the 
number of different detectors increased. For both raw response data and normalized data 

10 for the carbon black-polymer composite detectors, whose signal linearity vs. analyte 

concentration behavior permitted such an analysis to be performed in this test system, the 
same trend was observed for this particular array with this particular set of test analytes 
(as shown in Fig. 7B). 

At any array size and for any of the detector arrays investigated, subsets of 

1 5 the full array that can outperform the mean rf (n) performance of all possible 

combinations of elements having that specific array size can always be identified. 
Analysis provides an answer to the question, given that the task to best resolve on average 
all nineteen specific test vapors, whether the best performing array contains the full 
collection of available detectors or instead only contains subsets of each detector 

20 modality. The quantity that is evaluated is the average pairwise resolution factor, 

determined using the Fisher linear discriminant method, of the best-performing set of 
detectors in an array, rf max(n), at each number of detectors, n. Fig. 7C depicts the 
rf max(n) data for each of the three detector modalities using raw response data. For each 
of the three detector modalities, a larger number of detectors provides increased resolving 

25 power according to the rf max(n) criterion. 

In the case of the carbon black-polymer composite detectors, it is also 
possible to evaluate the concentration-independent performance using normalized data. 
As displayed in Fig. 7D, the optimum set of detectors obtained under these criteria once 
again is the full fourteen-detector array. 

30 Another goal might be the selection of the subset of detectors that 

maximized the smallest pairwise solvent resolution factor for the analytes in the test set. 
The data in Fig. 8A describe the ability of each array type to resolve the worst-resolved 
analyte pair as a function of the size of the detector array. For each of the three detector 
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types, the resolving power increases rapidly and then appears to plateau at large numbers 
of detectors. However, in none of the three cases, using raw response data to analytes at . 
fixed concentrations, does the addition of an extra detector diminish the resolving power 
of the array. 

5 The anay of eight tin oxide detectors maximizes the rf WO rst at 3.8 (for 

butylamine vs. triethylamine). This rfworst corresponds to approximately a 99% 
confidence in correctly distinguishing between the two analytes in a single presentation to 
the array. The array of twelve bulk conducting organic polymer detectors maximizes the 
rfworst at 1 .4 (for n-heptane vs. a,ct,a-trifluorotoluene). This corresponds to 

10 approximately an 84% confidence in correctly distinguishing between the two analytes on 
a single presentation to the array. The array of fourteen carbon black-polymer detectors 
maximizes rfworst at 23. For this best set of fourteen detectors, the worst resolved 
analyte pair is anisole and butyl amine. The rfworst of 23 essentially corresponds to a 
100% confidence in distinguishing between the two analytes In a single presentation to 

15 the detector array. 

Fig. 8B depicts rfworst va l ues f° r normalized data for the carbon black- 
polymer composite detectors. Again, the best resolution is maximized by the full 
fourteen detector array, although the data appears to plateau for sets containing more than 
eleven detectors. The full fourteen-detector set maximizes rfworst at 21 (for butylamine 

20 vs. n-pentane), which essentially corresponds to a 1 00% resolution probability. This full 
set of fourteen carbon black-polymer composite detectors contains polymers which vary 
widely in chemical properties, from the polar poly(N-vinylpyrrolidone) to the 
comparatively non-polar poly(ethylene-co-vinyl acetate), and includes halogenated and 
aromatic functional groups such as those in poly(vinylidene chloride-co-acrylonitrile) and 

25 poly(9-vinylcarbazole). The diversity in the detector component properties, which 
matches well with the diversity in the analytes, leads to the ability to resolve all the 
nineteen analytes in the test set quite effectively. 

The present invention also allows analysis in the situation where the task is 
known and can be specified precisely in advance of array construction. Members of the 

30 detector array can, in principle, be pre-selected to provide the optimum performance for 
the task of concern. The issue to be addressed is whether the optimal resolution of a 
specific solvent pair, at a fixed concentration, is obtained through use of the data 
produced by the full array, or whether it is advantageous, for that particular pair of 
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analytes, to use only the data produced by a subset of the full detector array. Ideally, all 
detectors would be orthogonal and one detector would probe exactly along the direction 
that is associated with maximizing the chemical differences between the two specific 
analytes. This single detector, and duplicates of it, should form the best "array" for 
5 resolving the two test analytes, under controlled conditions, at fixed concentrations. 
However, since in reality the vapor detectors are not generally mutually orthogonal, nor 
do they probe exactly the chemical differences between a specific analyte pair, it is 
probable that a few specific, partially correlated, detectors will combine to resolve a 
specific analyte pair better than any single detector in the array. 

10 As a representative example, Fig. 9 depicts the results for resolving an 

arbitrarily chosen solvent pair, ethyl acetate from tetrahydrofuran, as a function of the 
size of the detector array, using only raw response data. All possible combinations of n 
detectors in the array are investigated for the fourteen detector carbon black-polymer 
composite array, the eight detector tin oxide array, and the twelve detector bulk 

1 5 conducting organic polymer array. Only elements from a given detector array are 
allowed to form the desired detector array (i.e., carbon black composites, tin oxide 
detectors, and bulk conducting organic polymers are restricted to be in mutually separate 
arrays for this analysis). For each specific number of detectors in the array, only the best 
resolution factors that resulted from the optimal array for the task of separating ethyl 

20 acetate from tetrahydrofuran, rfbest(n), are plotted. 

In each of the three cases, the ability of the arrays to resolve ethyl acetate 
from tetrahydrofuran increases as additional detectors are included. From Fig. 9, it is 
apparent that beyond approximately the first six carbon black-polymer detectors (#8-13 
from the table in Fig. 6), additional detectors provide minimal increases in the resolution 

25 of ethyl acetate from tetrahydrofuran. Five of these six carbon black-polymer composite 
detectors contain polymers with one or both of the ether and ketone functionalities 
present in the ethyl acetate and tetrahydrofuran analytes, which explains chemically why 
the detectors perform well in this specific task. The other detectors in the array produce 
response data along vectors that do not lie as close to the vector that best separated these 

30 two solvents in detector space. 

A similar metric is evaluated for all other solvent pairs as well. For each 
array class, using raw response data, the full array maximizes the resolution of each 
analyte pair. The best resolution factor that is obtained between any vapor pair, using any 
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set of detectors of the same class, is obtained for distinguishing n-pentane from toluene 
using the carbon black-polymer composite detector array which has rf = 556 using 
unnonnalized data for this task. Similarly, using normalized data from the carbon black- 
polymer composite detector array, the best resolved analyte pair is dichloromethane from 

. 5 a,a,a-trifluorotoluene with rf = 256. The largest resolution factors obtained between any 
analyte pair, using raw response data, for the conducting polymer or tin oxide detectors 
are 77 (for ethanol vs. isopropanol) and 148 (for chloroform vs. toluene), respectively, 
again for unnormalized data. 

One of the findings to which the present invention leads is that array 

10 performance increases, in general, as the number of sensors increases. The larger number 
of sensors is beneficial for resolving, on average, a generalized set of test vapors which 
might not all be known in advance of the array design. A larger number of sensors should 
increase the probability that the dimensionality of odor space is fully spanned by the 
array. Stated differently, the array should be able to probe many chemical differences 

15 between analytes. On the other hand, with the analytes are fixed , a smaller number of 
selected sensors could optically resolve between solvent pairs. An optimal number of 
sensors can be determined for a given set of analytes. 

The present invention is also capable of quantitatively pointing to the 
superiority of certain types of arrays for a given set of analytes. As between the tested 

20 arrays of carbon black-polymer composite detectors, the tin oxide detectors, and 
conducting polymer detectors, it is clear that the carbon black-polymer composite 
detectors outperform the other sensor arrays. However, combining arrays of different 
detector of different modalities should lead to synergistic performance for the resolution 
of many analytes. For example, bulk organic conducting polymer detectors might be 

25 placed in an array of surface acoustic wave (SAW) devices, tin oxide detectors, dye- 
impregnated polymer films on fiber optic detectors, quartz crystal microbalances 
(QCMs), or carbon black-polymer composite chemiresistors, for example. Besides the 
resolution performance, choosing between these various types of sensor arrays might 
require consideration of more practical matter, such as the signal-to-noise ratio of a given 

30 type of a signal transduction mechanism, the ease of manufacturing the arrays and their 
costs, and the integration with processing electronics. 

Hence the use of a numerical measurement of the resolution factor of a 
detector array according to the present invention can enhance the quantitative 
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understanding of the size of the sensor array, whether adding new detectors should be 
added, or whether certain chemical properties are underrepresented in a specific array, for 
example. The present invention helps the evaluation of other factors that are important in 
the design of arrays of vapor detectors. 
5 It should be noted that the Fisher linear discriminant technique described 

herein is one analytical approach for resolving between analytes. Another related 
approach is principal component analysis in which the data are transformed such that they 
project along mutually orthogonal vectors in n-dimensional space so that the first 
principal component is the vector that captures the most variance in the data set and 

10 higher principal components capture progressively less variance in the data. This 

statistical approach preserves the inherent weighing of the original data set in terms of 
molecular properties being probed by the sensors, but allows for improved visualization 
of the differences between the various patterns as they cluster in principal component 
space. In either principal component space or original, i.e., sensor, coordinate space 

1 5 (used in the Fisher linear discriminant technique), a statistical measure of the differences 
between two clusters of patterns can be obtained by determining the distance between the 
centroids of the data arising form repeated exposures to a given vapor and dividing this 
distance by the sum of the standard deviations of the two different pattern clusters 
projected along the vector that connects these centroid points. This metric is independent 

20 of the coordinate system chosen to display the data and is instead an inherent statistical 
property of the data set. 

From standard vector analysis, the mean response vector, x a , of an n-sensor 
array to analyte a is given as the n-dimensional vector containing the mean autoscaled 
response of each sensors, A a j, to the ath analyte as components such that 

25 X a = (Aal, Aa2, ... Aan) 

The average separation, I d | , between the two analytes, a and b, in the Euclidean sensor 
response space is then equal to the magnitude of the difference between x a and Xb. The 
noise of the sensor responses is also important in quantifying the resolving power of the 
sensor array. Thus the standard deviations , a^d and ab,d, obtained from all the individual 
30 array responses to each of a and b along the vector d, are used to describe the average 
separation and ultimately to define the pairwise resolution factor as 
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Further details of the present invention are described in the previously 
referenced U.S. Provisional Application No. 60/086,876, filed May 27, 1998 and U.S. 
Provisional Application No. 60/088,804, filed June 9, 1998, which applications are 
hereby incorporated by reference for all purposes. 

While the above is a complete description of the preferred embodiments of 
the invention, various alternatives, modifications and equivalents may be used. It should 
be evident that the present invention is equally applicable by making appropriate 
modifications to the embodiments described above. Therefore, the above description 
should not be taken as limiting the scope of the invention which is defined by the 
appended claims. 



WO 99/61902 1 g PCT/US99/11636 

WHAT IS CLAIMED IS: 

1 . A method of resolving different odors, said method comprising 
disposing a plurality of d sensors in an array, each sensor having different 

responses to said different odors; 

exposing said d sensors to first and second odors; 

generating first and second sets of data from each of said sensors, each set 
corresponding to said first or second odor, each data capable of being represented by a 
vector in a d-dimensional space; 

determining an axis in said d-dimensional space, said axis having the property that 
projections of said data points onto said axis in said d-dimensional space have optimal 
separation on said axis; and 

resolving said first odor from said second odor by said separation. 

2. , The method of claim 1 wherein d is greater than four. 

3. The method of claim 2 wherein d is at least eight. 

4. The method of claim 1 further comprising 
quantifying said resolving step by a separation amount 

5. The method of claim 1 wherein in said disposing step, each sensor 
has an electrical response to said different odors. 

6. The method of claim 1 further comprising 
disposing said sensors on a single substrate. 

7. The method of claim 6 wherein said substrate comprises 
2 semiconductor material . 

1 8. The method of claim 6 wherein said sensors comprise carbon 

2 black-polymer composite sensors. 

1 9. The method of claim 6 wherein said sensors comprise tin oxide 

2 sensors. 

1 1 0. The method of claim 6 wherein said sensors comprise bulk organic 

2 conducting polymer sensors. 
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2 11. The method of claim 8 wherein said sensors sense odor by a 

2 change in resistance responsive to exposure to said odor. 

1 12. The method of claim 1 wherein said sensor disposing step 

2 comprises 

3 disposing a plurality of different types of sensors in an array. 

1 13. The method of claim 12 wherein said disposing step further 

2 comprises selecting at least one type of sensor from the group comprising surface 

3 acoustic wave devices, tin oxide detectors, conducting organic polymers, dye- 

4 impregnated polymer films on fiber optic detectors, polymer-coated micromirrors, quartz 

5 crystal microbalances, electrochemical gas detectors, chemically sensitive field-effect 

6 transistors, carbon black-polymer composite cherniresistors, micro-electro-mechanical 

7 system devices, and micro-opto-electro-mechanical system devices. 

1 14. The method of claim 13 wherein at least one type of sensor 

2 comprises carbon black-polymer composite sensors. 

1 15. The method of claim 1 3 wherein at least one type of sensor 

2 comprises tin oxide sensors. 

1 1 6. The method of claim 1 5 wherein at least one type of sensor 

2 comprises bulk organic conducting polymer sensors. 

1 17. The method of claim 1 6 wherein a second type of sensor comprises 

2 carbon black-polymer composite sensors. 

1 1 8. A system for resolving different odors, said system comprising 

2 an array having a plurality of d sensors, each sensor having different electrical 

3 responses to said different odors; 

4 an array interface receiving first and second sets of data from each of said sensors 

5 upon exposure of said d sensors to first and second odors, each set corresponding to said 

6 first or second odor; and 

7 a computer unit storing said first and second sets of data, said computer unit 

8 processing each datum as a vector in a d-dimensional space, determining an axis in said 

9 d-dimensional space, said axis having the property that projections of said data points 
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1 0 onto said axis in said d-dimensional space have optimal separation on said axis and 

1 1 determining said optimal separation to resolve said first odor from said second odor by 

12 said separation. 

1 1 9. The system of claim 1 8 wherein said computer unit has a unit upon 

2 which said optimal separation is displayed. 

1 20. The system of claim 1 8 wherein said computer unit has a unit upon 

2 which said optimal separation is stored. 

1 21 . The system of claim 1 8 wherein each sensor has an electrical 

2 response to said different odors. 

1 22. The system of claim 1 8 wherein said sensors are disposed on a 

2 single substrate. 

1 23. The system of claim 22 wherein said sensors comprise carbon 

2 black-polymer composite sensors. 

1 24. The system of claim 22 wherein said sensors comprise tin oxide 

2 sensors. 

1 25. The system of claim 22 wherein said sensors comprise bulk organic 

2 conducting polymer sensors. 

1 26. The system of claim 1 8 wherein said d sensors comprise a plurality 

2 of different types of sensors. 

1 27. The system of claim 26 wherein at least one type of sensor is 



2 selected from the group comprising surface acoustic wave devices, tin oxide detectors, 

3 conducting organic polymers, dye-impregnated polymer films on fiber optic detectors, 

4 polymer-coated micromirrors, quartz crystal microbalances, electrochemical gas 

5 detectors, chemically sensitive field-effect transistors, carbon black-polymer composite 

6 chemiresistors, micro-electro-mechanical system devices, and micro-opto-electro- 

7 mechanical system devices. 

1 28. The system of claim 26 wherein at least one type of sensor 

2 comprises carbon black-polymer composite sensors. 
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1 29. The system of claim 26 wherein at least one type of sensor 

2 comprises tin oxide sensors. 

1 30. The system of claim 26 wherein at least one type of sensor 

2 comprises bulk organic conducting polymer sensors. 

1 31. The system of claim 30 wherein a second type of sensor comprises 

2 carbon black-polymer composite sensors. 

1 32. A method of resolving different odors, said method comprising 

2 disposing a plurality of sensors for exposure to said odors, said sensors 

3 including at least two different types of sensors; 

4 exposing said plurality of sensors to first and second odors; and 

5 generating first and second sets of data from each of said sensors, each set 

6 corresponding to said first or second odor; and 

7 resolving said first odor from said second odor from an analysis of said 

8 data sets. 

1 33. The method of claim 30 wherein at least one type of sensor is 



2 selected from the group comprising surface acoustic wave devices, tin oxide detectors, 

3 conducting organic polymers, dye-impregnated polymer films on fiber optic detectors, 

4 polymer-coated micromirrors, quartz crystal microbalances, electrochemical gas 

5 detectors, chemically sensitive field-effect transistors, carbon black-polymer composite 

6 chemiresistors, micro-electro-mechanical system devices, and micro-opto-electro- 

7 mechanical system devices. 



1 34. The method of claim 32 wherein one type of sensor comprises 

2 carbon black-polymer composite chemiresistors and a second type of sensor comprises 

3 conducting organic polymers. 

1 35. A method of resolving different odors, said method comprising 

2 exposing an array of odor sensors to a plurality of odors, each odor sensor 

3 having an electrical response to at least one of said odors ; 

4 determining the maximum relative electrical response for each odor and 

5 sensor; 
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6 determining the average maximum relative electrical response for the 

7 determined maximum relative electrical responses for each odor; 

8 determining the standard deviation for the average maximum relative 

9 electrical response for the determined maximum relative electrical responses for each 

10 odor; and 

1 1 computing the optimal resolution factor by the Fisher linear discriminant 

1 2 for each pair of odors. 

1 36. The method of claim 35 wherein the optical resolution factor rf 

2 equals 

^ 



4 where w is a projection vector, a %v/ is ^ e standard deviation for the average maximum 

5 relative response for odor a along vector w, Cb, w is the standard deviation for the average 

6 maximum relative response for the odor b along vector w, and d w is the distance between 

7 the average maximum responses for odors a and b along vector w. 



1 37. The use of Fisher linear discriminant optimum resolution factors 

2 from a set of data points obtained from exposures of an array of odor sensors to a 

3 plurality of odors to determine an optimum number of odor detectors for the array to 

4 resolve the plurality of odors. 

1 38. The use of Fisher linear discriminant optimum resolution factors 

2 from a set of data points obtained from exposures of an array of different odor sensors to 

3 a plurality of odors to determine a subset of the odor sensors to resolve selected odors to a 

4 desired degree of certainty. 

1 39. The use of Fisher linear discriminant optimum resolution factors 

2 from sets of data points obtained from exposures of arrays of different odor sensors to a 

3 plurality of odors to determine which array best resolves the plurality of odors. 
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